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Abstract 

We study effects of possible tachyonic perturbations of dark energy on the CMB 
temperature anisotropy. Motivated by some models of phantom energy, we consider 
both Lorentz-invariant and Lorentz-violating dispersion relations for tachyonic pertur- 
bations. We show that in the Lorentz-violating case, the shape of the CMB anisotropy 
spectrum generated by the tachyonic perturbations is very different from that due to 
adiabatic scalar perturbations and, if sizeable, it would be straightforwardly distin- 
guished from the latter. The tachyonic contribution improves slightly the agreement 
between the theory and data; however, this improvement is not statistically signifi- 
cant, so our analysis results in limits on the time scale of the tachyonic instability. 
In the Lorentz-invariant case, tachyonic contribution is a rapidly decaying function of 
the multipole number I, so that the entire observed dipole may be generated without 
conflicting the data at higher multipoles. On the conservative side, our comparison 
with the data places limit on the absolute value of the (imaginary) tachyon mass in 
the Lorentz-invariant case. 

PACS: 
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1 Introduction 



Recently, a number of suggestions have been put forward for explaining the observed accel- 
erated expansion of the Universe. Among them there are the presence of the cosmological 
constant, modification of gravity at ultra-large scales, existence of new light fields (for re- 
views, see, e.g., Refs. [U |2j [3l HJ E]). In the latter case, dark energy can be characterized 
by equation of state p = wp, where the parameter w is different from —1 and generically 
depends on time. In a simple version, dark energy is due to a scalar field (quintessence), and 
the parameter w obeys the constraint w > — 1, while w = —1 for the cosmological constant. 
However, the value of w may be strongly negative, w < — 1; this is the case for phantom 
dark energy. Existing data do not exclude also a possibility [7J, [8] that at relatively large 
redshift the dark energy is of quintessence type w > — 1, and later it becomes phantom with 
w < -1. 

Phantom energy violates null energy condition, which is usually a signal for instabil- 
ities. As an example, the simplest model of scalar field with the wrong sign of kinetic 
term [9] suffers from the presence of ghost (negative energy state) at arbitrarily high spa- 
tial momenta. This implies catastrophic vacuum instability. However, in phantom models 
that break Lorentz-invariance, violation of the null energy condition at cosmological scales 
(related to the property that w < — 1 for spatially homogeneous phantom field) does not 
necessarily imply unacceptable instabilities at shorter scales. This suggests that Lorentz- 
violating phantom theories may be viable. Indeed, models of this sort have been recently 
constructed dTj [T21 [T3]. A property of one class of these models [121 d3] is that there is a 
tachyonic mode in the perturbation spectrum about the homogeneous phantom background. 
This mode occurs at sufficiently small spatial momenta only, so that the time scale of the 
tachyonic instability may be roughly comparable to the age of the Universe. This is not 
particularly dangerous. 

It is conceivable that the existence of tachyonic modes at low spatial momenta is a fairly 
generic property of a class of phantom models: the violation of the null energy condition may 
show up precisely in this way. Therefore, it is of interest to study observable consequences 
of such models. In this paper we consider one of these consequences, namely, the effect of 
tachyonic modes on the anisotropy of CMB temperature. We adopt the phenomenological 
approach, and instead of using results obtained within a concrete model, we parametrize the 
tachyonic instability by a few parameters. We will consider models with a Lorentz-violating 
dispersion relation and Lorentz-invariant one. 

In the Lorentz-violating case, we parametrize the dispersion relation as follows: 



where a and M are constant parameters, and p is the physical spatial momentum. Our 
convention is that positive values of u 2 correspond to exponential growth of perturbations, 



u 2 = a\p\(M-\p\), 
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while the usual oscillatory behaviour occurs at negative u 2 . Thus, the parameter M equals 
to the momentum below which the mode is tachyonic. For given a, the parameter M 
determines also the time scale of instability. The parametrization ([1]) is chosen in accord 
with Refs. [T2"l [T3"] where similar dispersion relation has been found in a concrete model of 
phantom energy. The analysis presented below can be straightforwardly generalized to other 
forms of dispersion relation. 

In the Lorentz-invariant case the dispersion relation has the form 

uj 2 = M 2 - p 2 (2) 

In this case too, the tachyon mass M equals to the spatial momentum below which the mode 
is unstable, and M _1 is the time scale of instability. From our perspective, the important 
difference between the dispersion relations and (|2J) is that in the latter case the "fre- 
quency" is nonzero at p = and monotonously decreases as |p| increases, whereas in the 
former, the "frequency" vanishes as p = and has a maximum at finite |p|. This will lead 
to qualitatively different shapes of the contributions to the CMB anisotropy spectra. It is 
worth noting in this regard that the Lorentz-invariant model may be viewed as a represen- 
tative of a class of theories with tachyonic perturbations: analogous results would hold for 
Lorentz- violating models with dispersion relations similar to (121) . 

In the cosmo logical context, the dispersion relations and (T5]) are written as follows: 

u 2 (t) = oc4t(m k 



oj 2 (t) = M 2 — 



a{t) \ a(t) 
k 2 



a 2 (t) ' 

respectively, where k is the time-independent conformal momentum and a(t) is the scale 
factor. In the expanding Universe, a mode of given k is first normal, and after the physical 
momentum gets reshifted down to k/a = M, it becomes unstable. 

For M > Hq, where Hq is the present value of the Hubble parameter, the tachyonic 
modes in both models start to grow exponentially at times preceding the present cosmological 
epoch. The growth of the tachyonic perturbations gives rise to the growth of the gravitational 
potential $ generated by these perturbations^ [2] 



$(f,x) = ^^y"rfk$(t,k)e tkx + h.c., 
k) = A(k) exp j J 



u(t') dt' , (3) 



1 Refs. [T31 [H] consider the Lorentz-violating case, but that analysis is straightforwardly repeated in the 
Lorentz-invariant situation. 
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where A(k) is the amplitude of primordial fluctuations in the mode with conformal momen- 
tum k at the time tk when this mode becomes unstable. We will discuss the range of the 
amplitudes A(k) later on. We stress that the gravitational potential ([3]) is generated by 
tachyonic perturbations rather than inhomogeneities in the ordinary matter. 

In this paper we calculate CMB multipoles generated by the gravitational potential (^j. 
In the Lorentz-violating case, the mechanism we discuss gives rise to the contribution to 
the CMB anisotropy spectrum which is quite different from the standard spectrum coming 
from adiabatic scalar perturbations generated, e.g., at inflationary stage (for the latter see, 
e.g., Refs. [T51 ESI ELT|)- As we will see below, the tachyonic contribution leads to potentially 
observable features in the CMB spectrum at relatively low multipoles. On the other hand, 
there are hints towards the existence of deviations in the observed spectrum [6[ [18] from 
the predictions based on flat (Harrison-Zeldovich) or almost flat spectrum of primordial 
adiabatic perturbations. So, one is tempted to employ tachyonic instabilities to explain 
these deviations. We add the contribution due to the tachyonic perturbations to the standard 
contribution of the adiabatic modes and compare the result with the observed spectrum. We 
show that the tachyonic contributions slightly improve the agreement between the theory 
and data, but this improvement is statistically insignificant. So, our analysis only leads to 
limits on the parameters of the tachyonic perturbations. 

In the Lorentz-invariant case, tachyonic perturbations contribute to the lowest multipoles 
only. An interesting possibility here is that they may generate the entire observed dipole 
without getting in conflict with measurements of higher multipoles. 

This paper is organized as follows. We begin with preliminaries on the cosmological 
model in section [21 and then discuss in detail the growth of the gravitational potential due 
to the tachyonic modes in section [3], first in the Lorentz-violating model and then in the 
Lorentz-invariant one. In section H] we calculate the contributions to the CMB multipoles, 
again distinguishing Lorentz-violating and Lorentz-invariant cases. We compare our results 
with the data in section [5[ and conclude in section [6j 

2 The cosmological model 

The background space-time we consider in this paper corresponds to the "almost" standard 
cosmological model with the only special feature being that the accelerated expansion of the 
Universe is driven by phantom energy instead of the cosmological constant. The background 
metric is that of the spatially flat expanding Universe, 

ds 2 = dt 2 - a 2 (t)dx 2 . 
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The scale factor a(t) is determined by the Friedmann equation, which can be written as 
follows: 

'a(ty 



a(t) 



Hi 



a(t) J \a(t 



(4) 



where H Q is the present value of the Hubble parameter, a(t ) = a = 1 is the present value 
of the scale factor, the dot denotes the derivative with respect to cosmic time t, fl p and 
w p refer to phantom energy. We assume for simplicity that w p = p p / p p is independent 
of time; we will see in what follows that the effects we discuss are largely independent of 
w p , so this assumption is not restrictive. The values we use in this paper are Q m = 0.27, 
Q p = 0.73. According to observational data [6] the parameter w p belongs to the interval 
-1.38 < Wp < -0.86. 

When calculating the CMB multipoles, we work with conformal time rj instead of cosmic 
time t, 



v(t) 



dt 



o 



a(t) 

Equation (j3J), written in terms of conformal time, has the form 

da 



H dr] 



3 Growth of perturbations 



The gravitational potential of a mode in the tachyonic regime is given by eq. (131). The 
exponent 

t 

N{t,k) = [ diu{i) (5) 



k 



determines the growth of the potential in both Lorentz-violating model ([T]) and Lorentz- 
invariant one t^j. Properties of this function are different in the two cases, however. 

3.1 Lorentz-violating model. 

In the Lorentz-violating model, the function (jSJ) can be written as an integral over the scale 
factor, 



N{a 1 k) = ^Ja—^v \ da === , (6) 

-no J VdWtt m + n„a~ 3w p 



«k 
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where we have introduced dimensionless wave number 

k 

V= M- 

A few comments are in order. First, for a given mode the tachyonic regime begins when 
the integrand in becomes real, i.e., at a = v. Since a < ao = 1, the maximum value 
of v for tachyonic modes is v = 1, that corresponds to the modes that are entering the 
tachyonic regime today. In fact, because of the Hubble friction, the exponential growth of 
the tachyonic mode starts not quite at the time when a — v, namely, it starts at the time tk 
when the tachyonic "frequency" becomes comparable to the Hubble parameter, 

u(t k )~H(t k ). (7) 

The lower limit of integration in (jSJ) is the value of the scale factor at that time, a k = a(tk)- 
Second, the integral in (jfJJ) is a dimensionless smooth function of its arguments {a, u}, 
both of which do not exceed 1. Therefore, this integral is not parametrically large or small. 
On the other hand, the amplitude A(k) is small (see below), so the effect of the tachyonic 
instability may be considerable only if N(a, k) is sufficiently large. This can happen if 

y/aM > H Q . (8) 

We assume in what follows that this inequality is indeed valid. 

Third, the integral in (jBJ) is saturated near its upper limit, and because of (jHJ) this inte- 
gral is practically insensitive to the lower limit of integration. Therefore, it is an excellent 
approximation to set a k = u, that is, approximate the time tj, by the time at which the 
dispersion relation becomes tachyonic. In fact, by solving eq. (jTj) numerically, we have found 
that a>k = v with precision of order 10~ 5 in a wide range of values of v (1 > v > 0.05). By 
setting cik = v in ([6]) one finds that the dependence of N on the parameters M and a factors 
out 



N(a, v) = y/a—M{a, v) 



o 



Af(a, v) = ^fv J da 



\Ja — 



v 



—3w„ 



Fourth, the function jV(a, v) is obviously a growing function of the scale factor a, as 
shown in the left panel of Fig. [TJ This is especially relevant since the gravitational potential 
([3]) depends on iV exponentially. We present in the right panel of Fig. [1] the dependence of 
exp[iV(a, v)\ on the scale factor at v — 0.2; the plot is given for the values of parameter^! 



2 The effect of the tachyonic instability on the CMB anisotropy is sizeable if y/aM /Ho ~ 100 (see section 
hence our choice here. 
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Figure 1: Left: The growth function J\T(a,v) as a function of v at the values of the scale factor 
a = 0.8, a = 0.9 and a = 1.0 (dashed, dotted and solid lines, respectively). Clearly, N has a 
maximum at v ~ 0.2 and grows considerably as the scale factor approaches its present value a = 1. 
Right: The growth factor exp[N(a, u)] as a function of the scale factor a at v = f max ~ 0.2. The 
function exp[iV(a, v)] becomes comparable to its present value only at a > 0.95. Calculations are 
done for w p = -1, ^/aM/H Q = 100. 

such that tJolM/Hq = 100. It is clear that the major effect of tachyonic modes occurs at 
late times. 

Finally, let us discuss the dependence of J\f(a, v) on the parameter w p . We show in Fig. [2] 
the function J\f(a, v) at a = 1 for different values of this parameter. Clearly, the dependence 
on w p is weak; in what follows we take w p = — 1 for definiteness. 

The growth function jV(a, v) as a function of v has a pronounced maximum at i/ max = 0.2. 
The existence of a maximum of J\f(a, v) is a fairly generic property independent of the 
particular form of the parametrization ([I]), provided that w(p) vanishes at p = 0. Indeed, 
the tachyonic regime begins at k/a(t) = M. This implies, in particular, that M{a, v) vanishes 
at v — 1. Now, under the assumption that u;(p = 0) = one observes that M(a, v) vanishes 
at v = as well. Hence, Af(a, v) has a maximum at some intermediate v. In other words, 
modes of high conformal momentum have not entered the tachyonic regime yet. On the 
other hand, A/"(a, v) is small also at low conformal momenta, since the tachyonic regime 
begins too early for the corresponding modes. At that time the scale factor increases too 
rapidly; the function J\f(a, v) does not have enough time to grow by the epoch when the 
physical momentum k/a(t) gets close to zero and the growth terminates. The maximum 
growth occurs for intermediate momenta. We note in passing that the situation is different 
in the Lorentz-invariant model ([2]), since in that case u(p = 0) = M ^ 0, and the maximum 
growth takes place at the lowest spatial momenta. 

Besides the growth function N, the gravitational potential ([3]) is determined by the 
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Figure 2: Function N(a, v) at a = 1 for different values of the equation of state parameter 
of phantom energy, w p = —1.0, w p = —1.17 and w p = —1.33 (solid, dotted and dashed lines, 
respectively). 

amplitude A(k) at the time the perturbations enter the tachyonic regime. So, we have to 
estimate this amplitude at k/a = M. We assume that the tachyonic perturbations are 
generated from vacuum fluctuations. Hence, the amplitude is the Gaussian random field 
with zero expectation value, (A(k)) = 0. This field is completely determined by its two- 
point correlation function, which we parametrize as follows: 



where 5X is the primordial perturbations of the phantom field. Let us clarify the form chosen 
for the overall factor. We recall that the gravitational potential is generated by the pertur- 
bations of the phantom field, which in turn is assumed to cause the accelerated cosmological 
expansion. This implies, in particular, that the energy density of the homogeneous phantom 
field is of the order of the present critical density 



where Mx is the mass scale characteristic of the phantom field and X estimates the value 
of this field. By comparing eq. fllOp to eq. (1111) we obtain an estimate for this value, 




(9) 



p p ~ M 2 PL Hl 



(10) 



On the other hand, this energy density can be estimated as 



p p ~ M\X\ 




X 



M PL H 
M x 



(12) 
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Let us assume that all dimensionful parameters in the phantom Lagrangian are of order of 
the mass parameter M entering the dispersion relation (pQ), so that Mx ~ M. Therefore, 



5p p ~ M X X5X ~ MM PL H SX. 
Now, we make use of the Poisson equation for the gravitational potential, 



5p p 



M PL 



and find that at physical momenta of order M the gravitational potential is of order 

H -5X 



M PL ■ M 



This gives rise to the overall factor in (Q. The dimensionless function f(k) in (Q parametrizes 
possible deviations from the flat spectrum; it is likely that this function is model dependent. 
However, we have seen that the growth factor exp[N(a,v)] is peaked at u — v max . This 
implies that the integrals over conformal momenta, determining CMB anisotropies, are sat- 
urated in a narrow region of k. Hence, if f(k) is sufficiently smooth, our final results are 
insensitive to its shape, so we can set f(k) = C where C is a constant. This constant may 
be somewhat different from 1, as it may contain, e.g., a power of H /M; it cannot, however, 
contain extra factors involving the Planck mass. We will see that our final results are de- 
termined by the interplay between the large growth factor exp(iV) and small factor cx Mp\ 
in A(k), so a possible deviation of C from unity is unimportant in the end. Furthermore, 
we will see that the shape of the tachyonic contribution to the CMB anisotropy spectrum 
(the position of the maximum and width) is almost independent of M; the parameter M, 
as well as 8X, determine the overall magnitude only. Therefore, the constant C can be set 
equal to 1 by redefinition of the parameter M. In view of these observations we write for 
the two-point correlator, without loss of generality, 

^w>=3S^( k - k ')- < i3 > 

It is this expression that will be used in what follows. Note that the amplitude A(k) is 
small, so the linearized treatment of the problem is legitimate even for exponentially growing 
gravitational potential. 

Let us now discuss possible range of the amplitudes 8X. The lowest value is obtained by 
assuming that the perturbations of the phantom field are in the vacuum state just before the 
tachyonic regime sets in. Since we are interested in momenta of order M at that time, their 
amplitude is estimated as 5X ~ M. The largest conceivable primordial perturbations are 
those generated at inflationary epoch: their amplitude is then comparable to the amplitudes 
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of other (nonconformal) light fields, 5X ~ H in fl/(2ii), where H in fl is the Hubble parameter 
some 60 e-foldings before the end of inflation [151 HSJ El- m view of the observational 
constraint Hi n ji < 1 • 10~ 5 Mpl, we estimate the maximum value of 5X as 5X ~ 10~ 5 Mp£. 
We will see that the interesting range of the parameter Mx ~ M is roughly M < 10 3 H , so 
the latter amplitude is small compared to the background value (1121) of the phantom field. 
In what follows we present the results in the two extreme cases, 

H 5X M 



MM p L M PL 
and 

lO" 5 ^- (15) 



MM p L M 

The first case is obtained for vacuum amplitudes of fluctuations 5X before they get tachy- 
onic, and we neglected the relatively mild difference between M and H . The second case 
corresponds to the generation of primordial phantom perturbations at inflation with the 
maximum possible Hi n fl. We will see that the fact that these values differ by many orders 
of magnitude, the properties of the contributions to the CMB spectrum are qualitatively 
similar in the two cases. 



3.2 Lorentz- invariant model. 

The analysis performed in section 13.11 is straightforwardly repeated in the case of Lorentz- 
invariant dispersion relation (T5]). There are two properties, however, that make this case 
different from the Lorentz-violating one. The first property has to do with the behaviour of 
the function J\f(a, v) that enters the growth factor, 

N(a,v) = ^-N{a,v). (16) 
In the Lorentz-invariant case one has 



AV,,) / dh F 1« ^ T= • (17) 
V a\/\l m + iLa 6w p 



Unlike in the Lorentz-violating model, J\f(a, v) as a function of v monotonously decreases 
as v increases, so the maximum of Af(a, v) is at v = 0. It is therefore instructive to find 
J\f(a, v) at small v. To this end, we write for the integral f[T7]) 
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where 

A/q = I da 



Mi = I da (19) 
V ayn m + Q p d~ 3w p 



M 2 = - I da - _ (20 ) 



The integrals f|T9|) . fl20l) are convergent, since at large a one has 

,2 



Va 2 — v 2 — a ~ . 

2a 

Furthermore, the latter expression shows that at small ^ the integral ( l20l) is of order v 2 . The 
integral (TT91 can be written as 



IX) 

M = z/ 3/2 j dy 



V - 1 - y 



^/y^/Q m + Q p u- 3w py- 3w P ' 



Since W7 P < 0, the term with Q p in the integrand can be neglected at small z/, and then 
the remaining integral is straightforwardly evaluated. Finally, the integral Ao is readily 



calculated ato w p = — 1 , 



K = I la - £ = * ( Arcsinh,/^ - Arcsinh,/^ 



2 ArcsinhWa 3 ^- ^=z/ 3 / 2 + 0(z/ 9 / 2 ) 



With all contributions included, we obtain finally that at small v 



Af(a, v) = ArcsinhJa^ - -±=v^ + 0(A (21) 



CO 



Note that the expressions (fT6l) and fT2~TT) imply that the tachyonic contribution is non- 
negligible only for M > Hq. This is the analog of the inequality (JH])- 



5 Like in the Lorentz-violating case, the dependence on w p is weak. 
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It is clear from eq. (l2Tj) that unlike in the Lorentz-violating case, the growth factor 
exp[N(a,u)} is peaked at v = 0. This is due to the fact that low-momentum modes are 
tachyonic already at early times, and their "frequency" u = M is not small. However, the 
modes of very low spatial momenta k = Mv contribute to the monopole CMB harmonic only. 
This contribution strongly depends on the primordial spectrum of the tachyonic perturba- 
tions, i.e., on the shape of the function f(k) in ([9]); this is the second special property of the 
Lorentz-invariant model. Now, the monopole contribution merely renormalizes the average 
CMB temperature; it is not directly measurable and will not be discussed in this paper. 
Multipoles with I ^ are less model-dependent: we will see in section H~2l that the integrals 
over momenta are saturated in a relatively narrow region where N(a, v) is not damped, i.e., 
in the region where 

'§f ■ 

We assume that f(k) does not change much in an interval Ak ~ k for momenta belonging 
to the region (|23|) . Then for calculating the CMB multipoles with I ^ 0, we can still use the 
spectrum ( 1T31) . The above discussion of the primordial amplitude 5X holds in the Lorentz- 
invariant well. So, we will again concentrate on the two extreme cases ( JT4l) and 

(ED. 



4 CMB multipoles 



The contribution to CMB anisotropy we are interested in is generated fairly recently, when 
the tachyon-induced gravitational potential becomes sizeable. Therefore, the only phe- 
nomenon responsible for this contribution is the integrated Sachs- Wolfe effect. It is clear 
from Fig. [T] that in the Lorentz-violating model this effect operates at the late cosmological 
epoch beginning at z ~ 0.05 (a ~ 0.95). A similar picture holds in the Lorentz-invariant 
case. 

Let us use the standard notation for the temperature anisotropy 

e(-> = , 

where n is the direction of observation and T is the average CMB temperature at present. 
The integrated Sachs- Wolfe effect then reads (see, e.g., Ref. [19] ) 

r)(b(n 

(24) 



6(n) = 2 Id, d ^ 



o 



drj 



=n(r/o-r7) 



where r] is the present time, and in view of the above discussion we set the lower limit 
of integration equal to zero instead of the time of last scattering. The integrand in (T5] 
exponentially grows toward the present epoch. 
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CMB anisotropy is characterized by the multipoles 

m=l 
m=— I 

where ai m are the coefficients of the decomposition of the anisotropy over spherical harmonics, 

a lm = J dn9(n)y (ra (n). (26) 

Making use of eqs. ([3]), ffl3l) and (124]) and performing the angular integration one obtains 
the following expression for the multipoles in terms of integrals over conformal momenta: 

C = m \ S K f —Af(v) (27) 
■kM 2 M 2 pl J v n ' y ' 

o 

The quantity Aj(i/) is the analog of the power spectrum. It is expressed in terms of the 
integral over conformal time, 

10 

A t (u) = J dTLU^a^expiNla^^^MMu^o-r)], (28) 
>ik 

where 



is the spherical Bessel function of the first kind. The integrals ([28]) and (|27|) have different 
properties in Lorentz-violating and Lorentz-invariant models. 

4.1 Lorentz-violating model. 

We plot in the left panel of Fig. [3] the tachyonic contribution to the angular spectrum of 
CMB temperature in the Lorentz-violating model (|14p with a = 1 ■ 10~ 3 and M/H = 9770, 
using the standard quantity 

D, = <1±V Q . 
2tt 

A very similar plot is obtained for the Lorentz-violating model ffl5l) with a = 1 ■ 10~ 3 , 
M/Hq = 726. This is a very general situation: the tachyonic contributions to the CMB 
spectrum in models (|14p and (|15|) at the same value of a are virtually indistinguishable if 
the value of M appropriately scaled down by about an order of magnitude. The reason for 
that will become clear later, see eqs. (1351) . (136]) . The analysis of this section applies to both 
models (fl4l and (TT5T) : we will illustrate our results using the model (TL4"]) for definiteness. 
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Figure 3: Left: The tachyonic contribution to the angular spectrum of CMB temperature in the 
Lorentz-violating model (|14p with a = 1 • 10 -3 and M/Hq = 9770. Right: Same, but for a = 1.0 
and M/Hq = 315. 

It is clear that the spectrum shown in Fig. [3] has a rather narrow maximum at I = Z max 
(I max ~ 7 in this example). We will see below that the position of the maximum is determined 
solely by the parameter a, and that Z max grows as a decreases, see eq. (1361) . 

For a — 1, the spectrum is shown in the right panel of Fig. [31 It is clear from this figure 
that at large enough a, the tachyonic perturbations contribute to the lowest multipoles only. 
It is also clear that the dipole and quadrupole components differ by one order of magnitude 
only. Hence, we disregard in what follows the dipole component, as for realistic quadrupole 
it is negligible compared to the observed dipole, which is supposedly due to the motion of 
the Earth in the CMB reference frame. We note in passing that the situation is different in 
the Lorentz-invariant model: it is meaningful to discuss the dipole component in that case, 
see sections 14.21 and 15.21 

The multipoles C\ can be calculated analytically in the regime a(M/H ) 2 ^> 1, see (J5J). 
As we discussed in section |3T| the function exp[iV(a, v)\ rapidly grows with the scale factor, 
and as a function of v it has a peak at z/ max ~ 0.2. Thus, the main contribution into the 
integral (|2"51) comes from late times, while the integral (|2~T|) is saturated at v w z/ max (it is 
important at this point that i/ max is different from both and 1). This means, in the first 
place, that the lower limit of integration in (1251) may be set equal to zero. Second, one can 
make use of expansion of the function N(a, v) near a = 1, 



N[a(ri), u] = Nil, v) - u(v) ■ ( Vo - V ), 



(29) 



and for v = v. 



max — 



0.2 we have 



max 



) = 0.39^ 



M 
H~o' 
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Third, one can set the pre-exponential factor u(rj)a(rj) in ([28]) equal to its value at a = 1. 
In this way we obtain 

uMr/o 



exp [N(l, u)] / dxexp —\fa 



v 



-x ji(x) 



(30) 



where we introduced the integration variable x = i>M(r] — if). Finally, let us recall that 
i] ~ 1/Hq. Then the upper limit of integration in (l3"Uj) is uMt/q ~ z/ max M/if 3> 1, so the 
integration may be extended to infinity Then the integral in (130]) is calculated by making 
use of the formula [201 



J l+ i(x) 

dx 2_ — eX p (—yx) 



[1 + y 



1 T T{1 + 1)P_! * 



7 



7' 



r(i + i) 



i + f»)ir(z + |) 



1 + 2 



11, 3 l-^ 

-,-,/ + -, 

2'2' 2' 2 



,(31) 



where P _[ 2 is the Legendre function, F is the hypergeometric function^, and 



7 



a(l-i/) 
^ + ail — v) 



Because of the exponential dependence on v of the factor exp[iV(l, v)\ in (130]) . the integral 
in (|27]) can be evaluated in the saddle-point approximation, and we obtain 



D, 



2H 2 5X 2 
M 2 Mp L v max 



7r|iV"(l,z/ max )| 



exp(2iV(l,z/ max )) •/(/ + !) 



rq + i) 
r(/ + 3/2) 



x 



1 - z 



2 \l + z 



1+1/2 



1 1 



F -,-,/ + -, 
2 2 



3 1-2 



where z = z(u max ), N"(l,u) is the second derivative with respect to v. At u n 
have 



T(l + 3/2) 
2v^ 



/(/ + !) f l-z 
Vl + 4« \l + 2 



i+l/2 



11, 3 1-2 

-,-:/ + -, 

2'2' 2' 2 



(32) 



0.2 we 



(33) 



A 



Vl +4a 
9.4 • « 3 / 4 



0.78^— 

-no 



(34) 



4 The second representation of the integral (|3ip is convenient because the hypergeometric function is only 
slightly different from unity at < z < 1, which corresponds to < a < oo. 
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This is the desired analytical expression for the multipoles; it works with the precision which 
is certainly sufficient for our purposes.. 

At a <C 1 and I > 1, the expression (1321) simplifies to 



2aH$5X 2 1 - z/ max / 1 , , NW/,-n ( or r I 1 ~ ^™*\ 

= M 2 Mp L ^ v ^"' 1 -^ ( ( } } r 2 '^v^rj • 

(35) 

Corrections to the latter formula are (9(1//) and they are numerically small even for I = 2. 

A few comments are in order. First, it is seen from (I3"3"|) that the dependence on / and on 
M has factorized. Therefore, the position of the maximum Z max and the width of the peak 
in the spectrum depend on the parameter a and do not depend on the parameters M and 
5X. In particular, for a C 1 the maximum of the function fl35|) is at 



^max = 7r~7=l/ i UUiaX 1 - 7~7= ~ L ( 36 ) 

The overall magnitude of the spectrum depends on a, M and 5A. The dependence on M is 
exponential, since N(l, z/ max ) oc M. This justifies the use of ( TTBl for the primordial spectrum 
of perturbations of the gravitational potential. 

Second, the exponential dependence of C\ on Z has the following interpretation. The 
problem has the characteristic time scale r(k) ~ cj" - 1 = (y/EMy/u(l - v))- 1 < Fq- 1 . This 
scale determines the time of the development of the tachyonic instability in a mode with 
momentum k. Since the growth function N(y) has a maximum, the relevant modes have 
momenta near /c max = u max M. Therefore, the gravitational potential is small at distances 
r 3> r(/c max ) along the light cone emanating from the observer (the tachyonic instability 
has not developed yet). On the other hand, at distances r < r(A; max ) the gravitational 
potential is almost constant in time (according to (jEJ), the expansion of the Universe has a 
negligible effect, while the tachyonic instability gives rise to mild growth of the potential). 
In other words, the gravitational potential at r < r(fc max ) is the superposition of random, 
time-independent waves with almost constant amplitude and almost constant wavelength 
2vr//c max . At r > r(/c max ) the amplitude of these waves decays as exp[— r/r(/c max )] as r 
increases. The period of a wave located at distance r is seen at an angle A6 r ~ 27r/(rA; max ). 
Hence, this wave contributes to the multipoles with I ~ T/c max . The multipoles ai m are not 
exponentially suppressed for r < r(/c max ), i.e., / < l/[fc max r(fc max )], and are exponentially 
small in the opposite case. Recalling (|25|) . one finds that this behaviour of ai m leads to the 
following dependence of Ci on I, 



n I 21 \ ( oi r 11 ~ v " 

Ci oc exp — = exp —2l\/a\ 



^max^~(^max) / \ V 

in complete agreement with (|36l) . 
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The fact that for relatively large a sizeable contributions are obtained by the lowest 
multipoles only (see Fig. [3]) can be seen directly from (|24|) . Indeed, inserting ([3]) into (1241) . 
making use of (|29|) . and integrating over time, we obtain 



S(n) ~ d d kF(k,u)A(k) == + h.c, (37) 

where F(k, u) is a smooth function independent of the direction of k. Recalling that v ~ 
^max = 0.2, we have y/aty/(l — v)jv > 1 for sufficiently large a. In that case the denominator 
in fl371) can be expanded in a series in (k ■ n), which just corresponds to the expansion in 
spherical harmonics. The l-th harmonic is thus suppressed as (y/ciyj (1 — v)/v)~ l . This is 
relatively mild suppression, in accord with the right panel of Fig. [31 At \fa\/ (1 — v)jv < 1 
the expansion of the denominator in (1371) is not legitimate, and one has to perform more 
sophisticated analysis leading to (1331 . 

4.2 Lorentz- invariant model. 

In the Lorentz-invariant case we make use of the expression ([2"Tj) for the function J\f(a, v) to 
calculate the integral (1281) . We find 



lK J 2r(/ + 3/2)V2/ P ^3^ 0v /H; V^m ^ov^ J 

When obtaining this expression we used the fact that v is small (see ( 123|) ). again extended 
the time integration to infinity in the same way as we have done after eq. (1301 and made use 
of in eq. (13T|) . We also kept the leading terms in v in the expression (138]) . 

Inserting (I3"8"j) into the integral in the expression (1271) for multipoles, and changing the 
integration variable, we arrive at the following integral: 



41/3 4«/3 



H^yVtm \ f 4//3-1 _r / Hoy/Q m \ ( 4/ CM 





where j(/3,x) is an incomplete T function. We notice that the second argument of this 
function 



CM 



x 



is large, x ^> 1, otherwise the overall factor in (1271) makes the effect we discuss negligibly 
small. Hence, for (3 = 41/3 < i we use the approximation j(P,x) = r(/5), and obtain the 
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final formula for the multipoles 



D, 



H 2 SX 2 



1(1 + 1) 



3nM 2 M 2 L 2 2 '" 1 



cxp 



4M 1 
3H ,/U 



Arcsinh 



r(/ + 3/2) 



2CM 



3 

r 



This formula is valid for relatively low multipoles, I < CM / (Ho\/Tl r 
the behaviour of the incomplete V function is 



(40) 

In the opposite case 



l(P,x) ^ -je x 

Hence, the multipoles at large I are negligibly small in the interesting case of large x. 

Because of the exponential factor in (|40p . the low multipoles may be fairly large at large 
enough Mj H . On the other hand, a large value of this parameter implies that the multipoles 
rapidly decay with /. Consider first the case (fT4l) of small primordial amplitude SX. As we 
discuss in section 15721 the interesting range is M/Hq ~ 100, in which case the multipole -D/+i 
is suppressed by a factor of about 10~ 3 as compared to D\. So, it is sufficient to consider the 
dipole and quadrupole anisotropies only. With our parametrization, these are, respectively, 



D 2 



6.4- 10~ 125 



4.1 ■ 10 



-126 



/ M N 


2 /3 




) exp J 




2/3 




) exp ( 



1.98 



1.98 



H~o) 

M 



15.7 



M\ 4/3 
H ) 



(41) 
(42) 
(43) 



These expressions are in agreement, within 10%, with the values that we obtained numeri- 
cally. 

Let us now turn to high primordial amplitude SX, case (TI5]) . In that case, the interesting 
range is M/H ~ 10, and 

14/3 



Do 



2.7- 10 



-14 



HoX 
M J 



exp 1.98 ■ 



M 



while eq. ( |43j) remains valid. 



5 Comparison with the data 

Overall, the data on the anisotropy of CMB temperature are in good agreement with the 
standard picture of adiabatic scalar perturbations whose primordial spectrum is close to the 

5 The reason for the 10% discrepancy is the omission of the term of order v 2 in the expression (I2ip for 
Af(a, u); note that 7V(a, v) enters the final result exponentially. 
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Harrison-Zeldovich one. Still, the observed angular spectrum may possibly show deviations 
whose nature is unclear. Our analysis was partially motivated by the desire to understand 
whether these deviations may be due to the contributions coming from the tachyonic pertur- 
bations of phantom energy. As we discuss in this section, the deviations cannot be explained 
in this way. So, our analysis enables us only to place limits on the parameters of the tachyonic 
perturbations. 

5.1 Lorentz- violating case. 

We begin with the Lorentz- violating model, and consider a wide range of the parameter a, 

2.5 ■ HT 7 <a< 1.0. (44) 

This range is representative: at a > 1 the main effect is in the lowest multipoles, whereas 
at a = 2.5 ■ 10~ 7 the contribution is peaked at I ~ Z max ~ 500, see (13"oT) . As we pointed out 
above, we do not consider the dipole anisotropy in the Lorentz-violating it is much 

smaller than the observed dipole anisotropy supposedly originating from the motion of the 
Earth. 

In our study we used the data on multipoles given in Ref. [21] and organized as a table 
"C; vs. Z." These are not combined in bins, unlike the data usually presented. 

In the model we study, the multipoles are the sums of two terms, one due to adiabatic 
scalar perturbations and another due to tachyonic modes, 

Q = C\ ad) + C® . (45) 

Making use of the code CMBFast [22j, we calculated the angular spectrum Cj generated 
by adiabatic scalar perturbations for various values of the spectral index n s in the range 
0.8 < n s < 1.5, in the standard cosmological model with the following values of parameters: 
the Hubble constant H = 72 km • s _1 ■ Mpc -1 , baryon plus CDM contribution to the 
present energy density fl m = 0.27, contribution of hot dark matter Qhdm = 0, dark energy 
contribution Q p = Q\ = 0.73, 4 He abundance Y = 0.24, number of massless neutrino species 
N v = 3. We assumed that the tensor perturbations are absent. The second term Cj in (1431) 
was calculated by making use of the analytical expression (1331) . 

To compare the model with the data, we analyzed the difference between the measured 
and calculated multipoles, 

The study of the moments and correlation properties of e\ has shown that they are inde- 
pendent and their average is zero within statistical error. Furthermore, an x 2 estimate has 
shown that with 95% probability the distribution of e; is Gaussian. 
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To obtain the parameters of the theoretical spectrum, we used maximum likelihood 
method with the likelihood function 

F(e\6) = J] f^\ h 

Z=2,600 

em 



f(ei\0) = exp 



2af 



Here 9 is the set of four parameters: the spectral index and amplitude of adiabatic pertur- 
bations, the amplitude Aq and parameter a of the tachyonic perturbations. In view of (1361) 
and (|44|) . we included in our analysis multipoles with I < 600 only. For each a from the 
range (I44p we obtained the best fit values of the three other parameters. The variations of 
individual multipoles have been calculated with the use of the error estimations given in the 
third column of the table in Ref. [2Tj . 

We show in Fig. H] the best fit value of the tachyonic contribution as a function of the 
parameter a. As a measure of this contribution we use the maximum in the anisotropy 
spectrum generated by the tachyonic perturbations, 



D max — max 



2tt 



j i 



(46) 



This maximum is at I = / max (as an example, Z max ~ 7 in the left panel of Fig. [3]). It is 
clear from Fig. H] that at a > 10 -4 the best fit value is equal to zero, whereas at a < 10 -4 
it is considerably different from zero. This means that the tachyonic contribution improves 
the agreement between the theory and data. It is worth noting that the addition of this 
contribution moves the best fit value of the spectral index up from n s = 0.96 obtained in 
Ref. [23]; in particular, for some values of a the best fit values of n s are larger than 1. 

This improvement, however, is not statistically significant. We show in Fig. [5] the maxi- 
mum likelihood function as a function of -D max at a = 1.8 ■ 10~ 6 . It is clear from Fig. [51 that 
even though the best fit value of -D max is nonzero, the difference of the likelihood function at 
the best fit value and at -D max = is small. The same is true for all values of a in the range 
considered. So, the data is consistent with the absence of the tachyonic contribution to the 
CMB temperature anisotropy. 

Thus, we can only place limits on the overall magnitude of the tachyonic contribution, 
.Ao? a t various values of a, which can then be translated into the limits on the physical 
parameter M/H . These limits, at 95 % confidence level, are shown in Fig. [6] for both cases 
of low primordial amplitude ffl4l and high amplitude ffl5|) . 



5.2 Lorentz-invariant case. 

As we have seen in section 14.21 in the model with Lorentz-invariant spectrum the CMB mul- 
tipoles generated by the tachyonic perturbations rapidly decrease as / increases. Therefore, 
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D n 




Q: 



Figure 4: The best fit value of the tachyonic contribution to CDM temperature anisotropy as a 
function of a. The parameter D max is defined in (|46p . The maximum of -D m ax is at a = 3.0 • 10 -6 , 
which corresponds to Z max — 143. 



only two multipoles — the dipole and quadrupole — are relevant for comparison with the 
data. The measured dipole component of CMB temperature is [21] d = 3.358±0.017 mK, and 
the direction in the Galactic polar coordinate frame is I = 263.86 ± 0.04°, b = 48.24 ± 0.10°. 
In the standard parametrization, the dipole anisotropy is 



exp 



1 

3^ 



m=l 



1.6 • 10" 



(47) 



while the quadrupole component is given by 



^2 



2 m=2 

-2 



2.9-10 



-ii 



(48) 



To obtain conservative limits on the parameter M of Lorentz-invariant tachyonic perturba- 
tions, we do not impose any priors on the contribution to the dipole anisotropy due to the 
motion of the Earth and on the quadrupole anisotropy generated by adiabatic perturbations. 
In the model with small primordial spectrum (j!4j) we make use of the expression (141 j) and 
take into account the fact that the dipole has 3 degrees of freedom. In this way we find that 
the analysis of the dipole anisotropy leads to the limit 



M 

— < 135.9 



at 95 % c.l. 



The limit coming from the quadrupole anisotropy is obtained by making use of eq. 
reads 



M 



< 136.4 at 95 % c.l. . 



(49) 
It 

(50) 
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F / F max 

0.8 - 



Figure 5: Likelihood function at the best fit values of the amplitude and spectral index of the 
adiabatic perturbations as a function of D max at a = 1.8 • 10~ 6 . The dashed line shows the limit 
on .Dmax &t 95 % confidence level. The dotted line corresponds to the best fit value of D max . The 
maximum value of the likelihood function is i*max = 

0.577. 

Interestingly, the limit coming from the dipole anisotropy is similar to that obtained from the 
quadrupole. One can turn this result around and speculate that the large observed dipole 
may be due to the tachyonic perturbations, with no contradiction to the data at higher 
multipoles. In our model the latter property is natural in the sense that the quadrupole and 
higher angular harmonics are small automatically. Another way to phrase this is to pretend 
that the observed dipole anisotropy is due to the tachyonic perturbations, i.e., equate (phtl) 
and ([17)1 . and then calculate the contribution to the quadrupole from (H2|) . This gives for the 
tachyonic contribution D 2 = 1.5 • 10 , which is safely below the observed value (148|) . The 
octupole is suppressed by another three orders of magnitude. So, our model would serve as 
an alternative to other explanations [25], [27J, [26] [28], [291 EH EI] of the large dipole component 
of the CMB anisotropy, if such an explanation were needed (for observational aspects of this 
issue see Refs. (32] [33] and references therein). 

In the model with large primordial perturbations (|T5|) . the situation is different. In that 
case, the strongest limit on Mj H is obtained from the quadrupole 

M 

— < 8.6 at 95 % c.l. . 
H 

As seen from eq. ( H31 . the dipole is not so much enhanced as compared to quadrupole; at 
M/H = 8.6 its value is Di ~ 1 ■ 10~ 8 , which is safely below the observed value (jUJ- On 
the other hand, the octupole and higher harmonics are still suppressed compared to the 
quadrupole by more than an order of magnitude. 
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Figure 6: Left: Upper limit on the amplitude Ao for a G [l.O -j- 2.5 • 10~ 7 ] at 95 % confidence 
level. Right: Upper limit on the parameter M of the tachyonic perturbations in units of the 
Hubble constant, at 95 % confidence level, in the model (fT3| (solid line) and model (fT5|) (dashed 
line). 

6 Discussion 

In this paper we have considered the effects on the anisotropy of CMB temperature due to 
possible tachyonic perturbations of dark energy. Because of the exponential growth, these 
perturbations may generate large gravitational potential $ at the recent cosmological epoch, 
and only at that epoch. This results in a sizeable Sachs- Wolfe effect. Note that the tachyonic 
perturbations we have discussed are unrelated to perturbations in baryons or dark matter, 
so their contribution to the CMB anisotropy does not correlate with the distribution of 
structure in the Universe. 

Our analysis was mostly motivated by the Lorentz-violating models of phantom energy. 
Hence, we have studied in detail the tachyonic perturbations with the dispersion relation ([T]). 
We have seen that their effect on the CMB angular spectrum has a pronounced maximum 
whose position depends on one of the parameters, a, and is practically insensitive to other 
parameters. It is expected that similar shape of the angular spectrum is characteristic to a 
wide class of models with Lorentz-violating tachyonic perturbations, as it is closely related 
to the fact that these perturbations become sizeable at late times only. 

We have also considered tachyonic perturbations with Lorentz-invariant dispersion rela- 
tion ([2]). In that case, an interesting possibility is that the largest contribution to the CMB 
ansotropy is received by the dipole component, and the angular spectrum rapidly decays 
with the increase of I. We have seen that even if the entire observed dipole anisotropy is 
attributed to the tachyonic perturbations, the quadrupole component generated by them 
is still consistent with the observational data. It is worth noting that this result should 



23 



be inherent not only in the Lorentz-invariant model, but also in Lorentz-violating models 
with tachyonic dispersion relations, provided that the "frequency" does not vanish at zero 
momentum and decreases as momentum increases. 

Our main conclusion is that even if perturbations of the tachyonic type exist in the Uni- 
verse, their contribution to the CMB anisotropy is small. Nevertheless, we do not exclude a 
possibility that growing precision of observations, and especially elaborate analysis of corre- 
lations between CMB anisotropy and structures in the Universe, may lead to hints toward 
the possible exotic property of dark energy, the tachyonic behaviour of its perturbations. 
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